Topological nature of polarization and charge pumping in ferroelectrics 
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Electric polarization in insulators is represented by the transferred charge through a shift of the 
Bloch wavefunctions induced by an adiabatic change of external parameters Q. It is found that 
this covalent/quantum contribution is determined nonlocally by the topological structure in the Q 
space. The condition for the charge pumping for a cyclic change of Q is also obtained. Applications 
of this picture to various organic ferroelectrics and BaTiOs are discussed. 
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Dielectric properties of insulators have been among the 
most important subjects in condensed-matter physics 0. 
For insulators, each electron is often regarded as "local- 
ized", and is treated as a point charge. Hence the clas- 
sical picture of polarization, i.e., the displacement times 
the charge, demands an ionic nature for ferroelectricity. 
On the other hand, it has recently been recognized that 
covalency of the electronic state also contributes to fer- 
roelectricity For this mechanism, extended Bloch 
wavefunctions in the crystal are essential. An interfer- 
ence pattern of Bloch wavefunctions is sensitive to a 
change of parameters such as atomic displacements, lead- 
ing to a shift of the electrons as a whole • 

Experimentally, several novel quasi-one-dimensional 
organic ferroelectrics have been found ^ |^ . In these 
ferroelectrics, the polarization is mostly perpendicular to 
displacements of atoms or molecules, and is much larger 
than that expected from the classical picture. In partic- 
ular, some of hydrogen-bonded organic co-crystals with 
neutral molecules show a large polarization arising from 
a slight change of the hydrogen position, which can never 
be understood in the classical picture |[5| . This motivates 
us to study one-dimensional models. 

In one dimension, the covalent contribution to the po- 
larization dP induced by an adiabatic ^change of external 
control parameters Q is expressed as 



(1) 
(2) 



dP = \e\A-dQ, 



dQc 



where A; is a wave number, and \4)i,{k)) and n^(k) repre- 
sent an eigenstate and a momentum distribution function 
for the band (including the spins) i/, respectively. Aa is 
a linear response of the covalent part of polarization, and 
a generalization of the Born charge tensor. The choice of 
Q is at our disposal, such as atomic displacements, elec- 
tric field, and pressure. Equation (j^J involves the Berry 
phase of the Bloch wavefunctions in the momentum 
space. Here one can see similarity to the quantum Hall 



effect and anomalous Hall effect [lo| . 

In contrast to the ionic contribution, it is difficult to 
have an intuitive picture for the electronic contribution 
arising from the covalency. Only the sophisticated first- 
principle band calculations have succeeded in estimating 
the polarization 0. Thus, it is highly desirable to ex- 
plore a quantum and intuitive picture of polarization in 
covalent ferroelectrics and to provide a guiding princi- 
ple for a design of such ferroelectrics, which is a main 
purpose of this Letter. The basic idea is to consider the 
space of external control parameters Q, instead of the 
momentum space which has been focused in the previ- 
ous works 0, 0, 0| . This reveals a hidden topological 
structure in the Q space, expressed in the form of the 
magnetostatics in this space. We show that a trajectory 
of band-crossing Q points constitutes a "circuit" with 
a quantized total "current". The vector field A{Q) in 
Eq. (0) then corresponds to the "magnetic field" gener- 
ated by the circuit. We can calculate a change AP in 
the polarization due to that of Q as a contour integral of 
AiQ). 

This framework is related to the quantized charge 
pumping discussed by Thouless in the following way. 
When the parameters Q are changed along a cycle C en- 
closing the "circuit" , an integer multiple of the charge 
e is transferred from one end to the other of the sam- 
ple, corresponding to the charge pumping. On the other 
hand, the usual polarization corresponds to a fraction of 
charge pumping where the C does not form a loop. This 
picture also facilitates to design insulators showing the 
quantum charge pumping. 

We take a one-dimensional two-band Hamiltonian; 

^ = E claHaa'ik,Q)Ck,.', (3) 

H{k,Q) = eo(fc,Q)+ J2 ho.{k,QW, (4) 

Q=l,2,3 

where Ck^a and c\ ^ represent the annihilation and cre- 
ation operators of an electron of a wave number k and an 
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internal degree of freedom u, which will be identified as 
orbitals later. CTq {a = 1, 2, 3) are the Pauli matrices. We 
have omitted the spin degrees of freedom, which can be 
taken into account by multiplying AP by two. Though it 
is a simplified model, it contains essential features of the 
nontrivial topology. Even in one spatial dimension, by 
regarding the parameters Qa as additional dimensions, 
the topological structure becomes equivalent to that of 
higher dimensional systems, and can be nontrivial. The 
eigenenergies are given by e±{k, Q) — eo{k, Q) ± h{k, Q) 
with h = \h\. Below we study geometrical and topologi- 
cal properties of A in insulating states where the chem- 
ical potential fi lies in the gap between the two bands; 
|eo(A:, Q) — n\ < h{k, Q) for any k. 
Equation (0) is then written in terms of h as 

A - - —h(—x —\ (5) 

J-7T/a \^ dk dQa J ' 

where h = h/h, and a is the lattice constant. Equation 
(O gives a geometrical meaning to the polarization; an 
infinitestimal change dP is given by a solid angle of a 
ribbon in the h space bounded by the two loops which are 
drawn by h{k, Q) and ft.(fc, Q -\- dQ) for —-K/a < k < n/a 
(Fig. nj . Therefore dP is maximized when a shift of the 
loop is "perpendicular" to the loop itself. 




FIG. 1: Change of polarization induced by an adiabatic 
change of Q is schematically shown in the Q space. 

Henceforth, we assume that Q has three components 
Q — (Qi: (92, Qs), with which we can use an analogy 
to magnetostatics; generalization to cases with many 
components is immediate. To investigate topological 
properties of A in the Q space, we introduce a "flux 
density" B = Vq x A from Eq. namely, Ba = 

--J-i%t^^0.M;ii A)- i^g^^^'^^i' ^(^) 

the following three properties, (i) B is nonzero only along 
strings corresponding to band crossings h{k, Q) = 0. 
Actually, it shows a (5-function singularity there, (ii) 



Vq • B{Q) = 0, and hence the string forms a loop or 
goes from infinity to infinity, as exemplified in Fig. |31 
(iii) The total flux of the string is quantized to be an 
integer. In the rest of this paragraph, we derive these 
properties. Since \h\ = 1, three vectors and 

are all perpendicular to h and hence are coplanar. 

This leads to B{Q) = 0. However, this does not apply 
for band-crossing Q points where there exists k satisfy- 
ing h{k, Q) = 0. Such band-crossing Q points form a 
curve, i.e., a "string" in the Q space. This leads to (i). 
(iii) is shown by considering the total flux of the strings 
surrounded by a loop C, 

/ dQ ■ A{Q) ^ f dSg- B{Q) (6) 

where S is the closed surface surrounded by the loop C. 
It corresponds to the wrapping number of the unit sphere 
with the periodic boundary condition both for k and Q 
directions. Namely, Eq. © is the Pontryagin index for 
the mapping from (k, Q) [Q is on the loop C) to h, and is 
always an integer, (ii) follows because A{Q) is a continu- 
ous physical observable with no gauge degrees of freedom 
and contains a singularity only on the strings. Practical 
calculations of B are explained in 12]. 

The property (iii) is related with the quantized charge 
pumping discussed by Thouless 0; when Q is changed 
along a closed cycle C within the insulating states, the 
amount of pumped charge is the integer © . It represents 
the linking number between C and the strings. 

Next we consider an asymptotic behavior of A{Q) near 
the string. Asymptotically, it is inversely proportional to 
the distance from the string, i.e., to the band gap. From 
Eq. (0), it can be shown that the direction of A{Q) near 
the string is the same as that of a "magnetic field" due 
to an "electric current" . Details will be presented in . 

It is useful to decompose A{Q) into transverse and lon- 
gitudinal parts, A{Q) = A^{Q)+A\Q), where Vq x^' = 
and Vq -yl* — 0. The relation of B{Q) to the transverse 
component A^{Q) is analogous to that of current to mag- 
netic field in the electromagnetism. Hence, A*{Q) can be 
constructed from B{Q) by the Biot-Savart law. In par- 
ticular, when the model has a form ft,(fc, Q) ~ f{k) + Q, 
the longitudinal component A'{Q) vanishes and A{Q) is 
completely described by the Biot-Savart law; namely, the 
dielectric response A of the system is determined by the 
string structure representing the band crossing. 

Let us show how this theory is applied in specific ex- 
amples. One simplest model contains s-like and p^j-like 
orbitals on a molecule, aligned along the x direction (see 
Fig. 121(a)). This is the simplest model for the hydrogen- 
bonded organic ferroelectrics Q. Let us introduce three 
transfer integrals tss, tpp, and tsp, and the Hamiltonian 
is given by 77 (fc, Q) = [{tss+tpp)/2] cos k +{tsp sin k)a2 + 
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{[{tss —tpp)/2] cosk)a3 + Q-a (see Fig. 121(a)). The corre- 
sponding string and the distribution of A{Q) is given in 
Fig. Ol(al). For example, Qi is the atomic displacement 
which breaks the inversion symmetry, such as the hy- 
drogen atom position in hydrogen-bonded ferroelectrics 
, or an external electric field along the chain direc- 
tion. The parameter represents an energy difference 
between s and Px orbitals: q'^^ = {cg — ep)/2. We can 
modify Qa by external means. For example, the pressure 
p can tune this energy difference, i.e., Q3 = Q'f^ + Upp 
(ttp! constant), which can be of the order of 10 meV with 
J) ~ 10 kbar |5|. The square of the external electric field 
Ez perpendicular to the chain also acts as Qa, because the 
s and Px orbitals have different polarizability. The string 
intercepts the Q3 axis at Qs — ±A, where 2A — \tss—tpp\ 
is a gap at fc = without the orbital energy difference. 
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FIG. 3: Geometrically different string structures of B (red) in 
the Q space for different classes of models in the text. In (al), 
we also show the vector field A. The loop C encompasses the 
flux string. The arrows on the flux strings show the directions 
of B in the strings. Total flux of each string is unity, except 
for the strings on the Q3 axis in (bl) and on the Q2 and the 
Qi axes in (b2), whose total flux is two. 



FIG. 2: One-dimensional two-band models for molecular or- 
bitals are schematically shown. 

As is evident from Fig.El(al), Qi, being perpendicular 
to the string, is most effective to produce the polariza- 
tion. The cases Qg"' > A and Q^^^ < A correspond 
to ionic and covalent states, respectively. From the Biot- 
Savart law, in the ionic state the polarization P decreases 
rapidly, as P ea{t'^Qi/{Q'f^)^) {t ~ tss,tpp,tsp, a: lat- 
tice constant) for Qg"' A since the contributions from 
the segments in the circuit cancels partially, whereas in 
the covalent state (Qg^"* < A) the polarization P is con- 
stant for small Qg"': P ~ ea{Qi/t). 

Let us discuss possibility of the charge pumping. As 
described above, the electric field and/or the pressure act 
as Qi = aiEx, Q3 = Qs"^ + asE'^ or Q3 = -I- app 
(ai, as, ap-. constants). Therefore, a cycle C in Fig. 13 
(al), namely a cyclic change of Ex and E^ (or p) leads 
to the charge pumping. Hence, to achieve this quantum 
charge pumping, the gap at k — 0, {Q^f^ ± A|, should be 
small enough to be comparable with the experimentally 
attainable values of a^E^ or Upp. Therefore it is promis- 
ing if we can design an insulating molecular crystal with 
a gap at fc = is ~ 10 meV. 

One can also consider a change of the transfer integrals 
due to distortion. This leads to the term (Qi cos A:)cri 
instead of Qitxi. With this change, the topology of the 
string also changes into Fig. |31(a2). 



Next, we consider the s- and p^-orbitals (Fig. [3 
(b)). The model Hamiltonian is H(k,Q) = [{tss + 



t 



pp)/2] cos/c 



Hits 



tpp)/2]cosk + A)o-3 + QiCTi 



((52sinfc)(T2 -l- Q3CT3. where Qi cx E^ mixes the s and 
Pz orbitals, Q2 represents the rotation angle of the Pz or- 
bital, inducing the transfer between the s and Pz orbitals, 
and Q3 oc E'^ changes the difference between these two 
orbitals. The string for this model is shown in Fig.|31(bl). 
Instead of Qiai term in the above Hamiltonian, one can 
introduce a transfer integral between s and Pz orbitals 
by breaking the symmetry between z and —z parts of 
the electron cloud. This leads to (Qicosfc)cri, and the 
resulting string is shown in Fig. 13 (b2). In both cases, 
the polarization response will be larger if the change of 
Q is along the direction around the strings. 

Even without two orbitals per one molecules, one- 
dimensional dimerized charge-ordered systems such as 
TTF-CA and (TMTTF)2PF6 gives the same string struc- 
ture as in Fig. ini(bl). The Hamiltonian reads 



n 



E [(V2)(cl< 



■i+l 



h.C.) 

+Qi{-iycic, + {Q2/2)i-incic.,, 



h.C. 



,(7) 



where Qi represents the charge density order while Q2 
the dimerization order (see Fig. 0)). This Hamiltonian 
corresponds to H{k,Q) = (tcosfc -I- (33)0-3 -I- QiCTi -I- 
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((52 sin fc)cr2 with a fictitious Q3, and the corresponding 
string is shown in Fig. 01 (bl). For TTF-CA, Eq. 
should be regarded as that for one of the spin and Qi is 
given by the energy difference of the HOMO of donor and 
LUMO of acceptor molecules. In TTF-CA, the molec- 
ular displacement changes the transfer Q2, giving rise 
to the polarization almost perpendicular to the displace- 
ment p|. In Ref. 0, the electron-phonon interaction, 
namely Q2, is estimated to be ^ 0.2eV while the site en- 
ergy difference is Qi ~ 0.4 eV. Because the cyclic change 
in the Q1-Q2 plane around the Q3 axis of Fig. |3| (bl) 
transfers the charge 4e, one can roughly estimate the po- 
larization P as P = (2e/7ra^)Arctan((52/Qi) — 0.3e/a^. 
This is consistent with the more elaborated calculation 
including the electron correlation (TMTTF)2PF6, 
on the other hand, is a quarter-filled system, and Eq. {T)) 
should be regarded as the spinless fermion model derived 
for the infinite U limit, describing the charge degrees of 
freedom. In this case, both Q2 and Qi originate from the 
spontaneous symmetry breaking, which occurs at a low 
temperature Therefore a sufficiently strong electric 
field along the chain could induce the rotation of {Qi,Q2) 
around the origin, which is similar to the charge density 
wave sliding. However, we note that the CDW has double 
periodicity and there is no phase degrees of freedom. 

(?+e2)/2 (f-e2)/2 



a 
-a 



FIG. 4: An ionic or charge-ordering dimer model for TTF-CA 
or TMTTF 



Ferroelectricity of perovskite oxides like BaTiOs in the 
tetragonal phase 0| can be essentially described in terms 
of a one-dimensional d-p model 0|, which is identical 
to Eq. JT)). The ionic contribution from atomic dis- 
placements, which is not included in our tight-binding 
Hamiltonian, is estimated as follows: the ionic config- 
uration of Ba2+, Ti2 *^9+ and O^-^^- together with 
each displacement 0.06 A, 0.12 A and —0.036 A along 
the chain ^3 gives a polarization of Pjon — 0.13e/a^. 
On the other hand, the ionicity leads to a rough es- 
timate of Q\l \/Qi + ^ 2.89/4, while according to 
Ref. [ill, the 2.4% change of Ti-0 bond lengths yields 
an 8.7% alternation of d-p transfers, i.e., Q2 ^ 0.087t. 
Therefore the covalent contribution is estimated to be 
Pcov — (4e/7ra^)Arctan((52/Qi) ^ O.lle/a^, which is 
comparable to Pion- Then, the total polarization is 
roughly estimated as P = Pion + Pcov — 0.24e/a^, which 
explains the experimental value (0.20-0.25)e/a^ 

The present argument can be generalized to higher di- 
mensions. In three dimensions, for example, by changing 



ky and , the string will form a three-dimensional man- 
ifold. In the simplest case the string will form a bundle, 
and A{Q) outside this bundle represents the polarization. 
More generic cases are to be discussed elsewhere |20j . 

In conclusion, we explored the topological and nonlo- 
cal nature of polarization in the space of the external 
parameters Q. The Q points where a band crossing oc- 
curs play a crucial role in determining the polarization 
of the insulating state. This clarifies also the condition 
for the quantized charge pumping, and will be useful for 
designing materials with large dielectric response. 
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